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Within the context of nonrelativistic potential models, we obtain several formulas 
(with varying degrees of rigor) relating the wave functions at the origin of the cc, be and 
bb S-wave quarkonium systems. One of our main results is a model-independent relation 
which seems to hold to within 3% for any reasonable choice of interquark potential and 
any choice of radial quantum number - - namely, |\E r 6c(0)| 2 — |\I/cc(0)| 1 ' 3 |\I/ 6 ;;(0)| ' 7 (the 
exponents are motivated in the text). One of the physical consequences of this result is 
the following relationship between heavy meson masses which we expect to hold at about 
the 10% level: Mb* - M Bc ~ (0.7)(M J/i} - M Vc )°- 65 (M r - M^J 0,35 . 
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Recent major advances in our understanding of the nonrelativistic limit of Quantum 
Chromo dynamics (QCD) have generated renewed interest in the calculation of the produc- 
tion and decay rates of heavy quark bound states IJ. However, these computations still 
contain numerous nonperturbative "parameters" which cannot as yet be accurately deter- 
mined analytically from first principles. They must either be fit to experiment, determined 
numerically from lattice simulations of QCD, or extracted by more phenomenolgical con- 
siderations (for example, from potential models) . One important class of such parameters 
is the wave function at the origin (WFO), Vl/(0), for an S-wave bound state of a heavy 
quark and anti-quark. (More generally, for a bound state with angular momentum £, we 
should consider the quantity d ty/dr e evaluated at the origin. However, in what follows 
we will concentrate on the case I = 0.) The WFO enters not only into the production and 
decay amplitudes for heavy quarkonium systems, but also into the determination of the 
hyperfine splitting in their mass spectra. Large QCD and relativistic corrections to the 
simple first order formulas which relate the WFO's to the above observables (especially for 
charmonium) make it difficult to extract precise information about the WFO's from the 
experimental data. Moreover, while there has been much recent progress (see, for example, 
0), lattice simulations of QCD are not yet accurate enough to be very useful. This leaves 
us at the mercy of a more model-dependent approach. 

In nonrelativistic potential model descriptions of heavy quarkonia, it is a simple nu- 
merical exercise to extract highly accurate values of the WFO's for any given choice of the 
static potential between a heavy quark and anti-quark. The problem lies in which potential 
to choose. Two potentials which yield very similar spectra for the heavy mesons can give 
very different WFO's for any given state. The WFO's seem to be quite sensitive to the 
global details of the potential, while the energy levels are only sensitive to the shape of 
the potential in the vicinity of the RMS radii of the states being studied. In other words, 
when this sensitivity of the WFO's is coupled with our ignorance regarding the details of 
the heavy quark potential, we again arrive at the sobering conclusion that an accurate 
determination of the WFO's seems beyond our reach. 

So what's a theorist to do? Well, the situation isn't quite so dire as we have made it 
sound. There are certain rigorous, qualitative statements that one can make about WFO's 
within the context of potential models which hold for very large classes of potentials - 
including those believed to be relevant for heavy quark systems. For example, in a concave 
downward potential V(r) (defined by V'(r) > and V"(r) < for all values of the 
interquark distance r), it can be shown that the square of the WFO for the IS state is 
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larger than that of the 2S state 0. (Most likely it is more generally true that the square 
of the WFO decreases monotonically with increasing radial quantum number, but this has 
not, to our knowledge, been shown yet.) Moreover, it can be proven from lattice QCD that 
the static potential between two heavy color sources is concave downward |4[] . An accurate 
enough determination of the WFO's for low-lying heavy mesons can be made from the 
experimental data in order to test this qualitative result, despite the size of QCD and 
relativistic corrections. And indeed it holds true in both the cc and bb systems. Another 
such result is that the square of the WFO for the IS state increases with the two-body 
reduced mass \i faster than linearly for a concave downward potential f|. (This result is 
also conjectured, but not proven, for higher states.) A comparison of the WFO's extracted 
from the cc and bb data is consistent with this theorem as well. 

Should we have panicked if either of these two results were violated by the WFO's 
obtained from the data? No, not necessarily. We could have attributed the discrepancy 
to (at least) two possible sources. First, the uncertainty in any WFO extracted from 
the data can be estimated from an educated guess at the size of the higher order QCD 
and relativistic corrections not included in the determination (as well as the experimental 
uncertainty in the measurement). It is possible that these corrections are larger than we 
expect, and therefore a WFO obtained with the truncated series is not as accurate as we 
thought. Second, it is possible that our naive picture of, say, the J/ip as a simple bound 
state of a c quark and a c quark in a relative S-wave interacting via a static potential is 
incorrect. For example, one may have to consider a dynamical treatment of excited glue 
inside the meson, or allow mixing with other angular momentum states and/or continuum 
states. There is a nice example of this latter possibility [Q. It is strongly believed that 
in a concave downward potential, the energy splitting between the (n + 2)S and (n + 1)S 
states is always less than the splitting between the (n + 1)S and nS states, for any n > 1. 
However, in the charmonium system, though the measured 3S-2S splitting is less than the 
2S-1S splitting, the 4S-3S splitting shows an increase over the 3S-2S difference. Accepting 
the truth of the above conjecture concerning energy splittings in a concave downward 
potential, how do we explain the experimental numbers? The answer is that the threshold 
for open charm production occurs between the 2S and 3S levels and induces substantial 
mixing of the 3S and 4S cc states with continuum states. Thus, what we experimentally 
identify as the 3S and 4S levels of charmonium actually have substantial DD and D*D* 
components, among others. A similar state of affairs occurs in the bottomonium system. 
Here, the measured 5S-4S splitting is greater than the 4S-3S splitting. The open bottom 
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threshold occurs between the 3S and 4S levels and causes substantial mixing of the 4S and 
5S bb states with continuum states. It is this mixing in both the cc and bb systems that 
seems to be responsible for the apparent violations of the above energy splitting conjecture. 
This interpretation is supported by a coupled-channel analysis in the charmonium system 
0. The only purpose in showing this example is to remind the reader that no matter 
how generally a certain result may apply within the context of nonrelativistic potential 
models, there are still assumptions that must be made in order to relate these potential 
model results to real observations. And these assumptions may not hold for all states in 
all systems. 

With this disclaimer behind us, we can now describe the results of this note. The 
starting point for our investigation is a recent paper by Eichten and Quigg || which tabu- 
lates the WFO's for various quarkonium states in an assortment of "successful" potential 
models. We list them below (in natural units). 

(1) The Martin potential |]: V(r) = Ar - 1 + C, where A = 6.898 GeV 1 ' 1 , m c = 1.8 GeV 
and nib = 5.174 GeV. 



(2) The log potential 0: V(r) = A£n(r/r ), where A = 0.733 GeV, m c = 1.5 GeV and 
m b = 4.906 GeV. 

(3) The Cornell potential Q: V(r) = -A/r + Br + C, where A = 0.52, B = 
(1/2. 34) 2 GeV 2 , m c = 1.84 GeV and m b = 5.17 GeV. 



(4) The Buchmiiller-Tye potential [|12|: This potential has a rather complicated position 



space form. It is linear at large distances and quasi- Coulombic at short distances. The 
deviations from pure Coulombic behavior reproduce the running of the strong coupling 
constant to next-to-leading order in QCD. The global shape of the potential is essentially 
determined by two parameters — namely, the QCD scale (in the modified minimal sub- 
traction scheme) A.jjg which the authors of [|12| fit to be 509 MeV, and the QCD string 
tension which they take to be 0.153 GeV 2 (motivated by the light meson data). The po- 
tential also depends on the number of "light" flavors rif. The authors take rif = 3 for 
r > 0.01 fm, and rif = 4 for r < 0.01 fm. The quark masses used are m c = 1.48 GeV and 
m b = 4.88 GeV. 

(The parameters C in (1) and (3) and tq in (2) are irrelevant for 1^(0) | 2 .) Eichten and 
Quigg treat the cc, bb and be systems. 
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The first thing that catches one's eye in glancing at these tables is the apparent 
randomness of the entries. Of course one can spot the aforementioned general trends - 
namely, for a fixed quark content, the square of the S-wave WFO decreases with increasing 
radial excitation, and for fixed quantum numbers the square of the WFO gets bigger as 
one goes from the cc to the be to the bb system (increasing reduced mass). However, 
besides these qualitative behaviors, no additional regularity is apparent. For example, the 
square of the WFO for the T changes by about a factor of 3 between the various potentials 
— potentials which yield basically the same low-lying spectrum! Things like the ratio of 
the i/j(2S) and the J/ij) WFO's, or the ratio of the T and the J/ip WFO's, also cover a 
large range of values. Can any additional statements about these numbers be made which 
possess some degree of model-independence? 

Before we address this question, we would first like to present our own version of the 
S-wave portion of Tables I-III in ||, which corrects some small numerical errors made 
there. For instance, it is well known that for power-law potentials V(r) = Ar a + C, the 
square of the S-wave WFO scales with reduced mass \i as ^u 3 /( 2 + a ) ||. This result can also 
be used for the log potential by putting a = 0. The results of || show a mild violation of 
this scaling (on the order of a few percent for all radial quantum numbers) which cannot be 
accounted for by rounding errors. Upon our redoing of the computations using the Runge- 
Kutta method for solving the nonrelativistic Schrodinger equation, we found results which 
satisfied the scaling laws (within rounding errors) for the log and Martin potentials, and 
typically disgreed with the results of || in the second significant figure. We also tested 
our program on potentials with analytically known WFO's, such as the Coulomb, linear 
and harmonic oscillator potentials, and obtained agreement with the exact results to at 
least six significant figures. We then ran our program on the other potentials treated in 
H, the Cornell and Buchmuller-Tye potentials, and found similar disagreements to those 
encountered in the log and Martin cases. It should be stressed however that these mild 
errors in no way affect the conclusions of ||. We just want numerical results which are as 
accurate as possible in order to test some approximate formulas relating different WFO's 
that we will derive later. 

For all of these potentials we display results for the ground state as well as the next 
five radial excitations. This goes a little further than the results in Many of these 
states lie above the threshold for open flavor production, and hence in a region where the 
WFO's have limited usefulness because of mixing with continuum states. However, these 
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numbers are still quite useful in checking the general validity of the analytic formulas which 
are to come. 

We have also added one additional potential to the table. 

(5) The Lichtenberg- Wills potential flT3|: V(r) = 8tt(1 - Ar) 2 /[(33 - 2n f )r£n(Ar)], where 
we choose A = 0.7 GeV, m c = 1.84 GeV and m& = 5.17 GeV. At short distances, the 
running of the strong coupling constant to leading order in QCD is reproduced if one 
identifies A = e"> Aqcd, where 7 is Euler's constant. We have also taken the number of 
light flavors nj to be three in all systems studied with this potential. 

Note that we have chosen the b and c constituent quark masses to be the same as for the 
Cornell potential. The parameter A was then chosen so as to obtain a low-lying meson 
spectrum reasonably close to that obtained from the Cornell potential. The differences in 
the WFO's between the Cornell and Lichtenberg- Wills examples are then basically due to 
the different shapes of the potentials outside of the region between the RMS radii of the cc 
and bb systems. A substantial difference can still be seen between the two sets of WFO's, 
again emphasising their sensitivity to global features of the interquark potential. 

We first became interested in finding regularities in these numbers after a comment 
made to one of us by Ira Rothstein. He was able to prove that, in nonrelativistic QCD 
in the limit as m& — > m c , the square of the WFO for any state in the be system is equal 
to the average of the squares of the WFO's of the corresponding cc and bb states, plus a 
correction of order S 2 where 6 = (mi, — m c ) / (mj, + m c ) [ |14|| . That is, 

|<M0)| 2 = (|* cc -(0)| 2 + |^(0)| 2 )/2 + 0(S 2 ). (1) 

We have shown that this is also true in an arbitrary nonrelativistic potential model. Indeed, 
one can prove a slightly stronger result: 

|<M0)| = (|^ cc -(0)| + |^(0)|)/2 + O(5 2 ). (2) 

A similar formula also holds for the geometric mean instead of the arithmetic mean: 

|^ c -(0)| 2 = |vl> cc -(0)||vl> 6 5(0)| + O(5 2 ). (3) 

These last two results can be easily demonstrated from perturbation theory in 5. However, 
even though these results are independent of the nature of the interquark forces, they are 
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unfortunately not very useful in real applications since the quantity 5 is approximately 1/2 
for reasonable values of mj and m c . The order 5 2 corrections in the above equations are 
therefore large, which a simple check using the numbers in Table I will show. 

What we want is a relation with the model-independence of Eqs.(l)-(3), but with 
much more quantitative accuracy. For the class of power-law potentials V(r) = Ar a + C 
discussed earlier, there is a very simple, exact relationship between the WFO's of the cc, 
be and bb systems. In order to derive this relation, we first recall that simple scaling 
arguments for the above power-law potentials tell us that for any reduced mass \i we 
have |* M (0)| 2 = f(n, a)(A^f^ 2+a \ where f(n, a) is only a function of the radial quantum 
number n and the power a. Using this fact alone, it is straightforward to obtain, for 
reduced masses Hi < ^2 < A*3 an d any fixed n, 

l*M a (0)| 2 = |*Mi(0)l 2(1 " ,) l*M3(0)l 2fl . (4) 

where q = tn{\ii j \x\) / £n(fis / Hi) . Choosing ^1 = m c /2, /i 2 = mi,m c /(m6 + m c ) and 
A 4 3 = m &/2, this becomes 

|* 6a (0)| 2 = l^eaCO)]^ 1 -^!^^)! 2 ^ (5) 

where q = £n(2mi> / (nib + m c ))/£n(mb/m c ). This result is nice not only for its simplicity, 
but also because it does not depend on any of the parameters (A, a and C) appearing 
in the potential. It depends only on the constituent quark masses nib and m c . It is easy 
to check this result on the log (a = 0) and Martin (a = 0.1) potentials in Table I. Since 
this formula has no dependence on parameters in the potential, we can also check it on 
the other examples in Table I. Of course it will no longer be exact in these cases since the 
above scaling law for ^^(O) | 2 is true (for all /u) only for power-law potentials. And these 
other potentials are far from being power-like. They each have a (quasi-) Coulombic nature 
at small r, motivated from one gluon exchange, and a (quasi-) linear behavior at large r, 
motivated by a stringy picture of confinement. In this sense, they are more "realistic" 
than the power-law potentials. In the intermediate r range containing the RMS radii 
of the heavy quarkonium states, they are quasi-logarithmic, just like the log and Martin 
potentials. But, though Eq.(5) is not exact here, we can still ask if it is a reasonably 
accurate approximation. 

The answer is yes. For every choice of n in Table I, the relation in Eq.(5) holds to 
within 4% (except for the IS state of the Cornell potential where it is off by about 7% ). 
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The least accurate results are obtained for the ground state. As n increases, the results 
get better. This is a substantial improvement over the accuracy of Eqs.(l)-(3). The only 
price that we have had to pay is the introduction of the constituent quark masses into the 
relation. It is interesting to note that the left-hand side of Eq.(5) is less than or equal to 
the right-hand side for each potential considered and each choice of n. Is it possible that 
this is always the case — at least for a wide class of potentials? A numerical study of 
numerous examples, as well as an analysis of the question within the context of various 
approximation schemes, has led us to the following conjecture: 

Conjecture: Consider a potential V(r) such that p(r) = rV"(r)/V'(r) is monotonically 
increasing with increasing r. Then for each choice of radial quantum number, and for 
reduced masses \i\ < ^2 < A*3, we have 

I^m 2 (0)| 2 < i^ 1 (0)| 2 ( 1 -^|^ 3 (0)| 2 ^, (6) 

where q is as in Eq.(4). For p(r) monotonically decreasing with increasing r, the inequality 
in Eq.(6) is reversed. 

Of course when p(r) is independent of r, V(r) is a power-law potential and the inequality 
in Eq.(6) is replaced by the equality of Eq.(4). One can think of 1 +p(r) as the "effective 
power" of V(r) at quark separation r. We will call a potential power increasing, or PI, when 
p(r) is monotonically increasing, and power decreasing, or PD, when p(r) is monotonically 
decreasing. Each of the non-power-law quarkonium potentials in Table I is PI (we have 
checked this numerically for the Buchmiiller-Tye potential), and satisfies the inequality in 
Eq.(6) with the appropriate choices of ^2 and ^3 — namely 

l<Mo)| 2 < \y c - c (o)\ 2{1 - q) \v b - b (o)\ 2q , (7) 

where q is as in Eq.(5). Indeed, all of the popular potentials used in quarkonium studies 
seem to be PI. But, unlike the concave downward property, we know of no QCD-motivated 
reason why this must be so. But we conjecture that Eq.(7) holds in the nonrelativistic 
limit of QCD, and in all realistic potential models. 

Actually, there is an even better result which is completely parameter-independent. 
To obtain this, we first note that q — £n(2mb/(mb + m c ))/£n(mi,/m c ) lies between about 
0.36 and 0.38 for any reasonable choices of m& and m c . However, as noted above, when 
substituted into Eq.(5) this yields WFO's for the be system which are too high. In the 
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context of the general form of Eq.(5), the potentials of interest seem to favor a slightly 
lower value of q. We have found that if q is simply taken to be 0.35 independent of the 
interquark potential and quark masses being considered, very accurate results are obtained. 
That is, we have 

l^o)! 2 ^!^)! 1 - 3 !^^)! - 7 . (8) 

Although no longer exact for power-law potentials, this simple formula holds to within 
2.5% for all cases in Table I. This is quite remarkable given the range of radial quantum 
numbers covered and the global differences in the potentials treated. We fully expect it 
to have a similar accuracy for any reasonable quarkonium potential. Though not on the 
same rigorous footing as the two qualitative theorems discussed earlier, it is reasonably 
well motivated by Eqs.(5) and (7) above. Moreover, it gives us a better quantitative 
understanding of the jumble of numbers in Table I. 

Can we extract any simple physical consequences of this result? Certainly it implies 
relationships between the production (and decay) amplitudes for the J/i/j, B c and T sys- 
tems. However, it is perhaps simpler to discuss the implications for the hyperfine mass 
splittings in these systems. To leading order in a s and f 2 /c 2 , the mass splitting (AM)- 
(for any fixed n) between the vector and pseudoscalar mesons composed of a quark of 
flavor i and an antiquark of flavor j (of mass rri\ and mj , respectively) is given by 

(AM)y = 327ra s (2^ iJ )|* iJ (0)| 2 /9m i m j , (9) 

where = m\m^/ '(mj + mj), and we have assumed the standard Breit- Fermi hyperfine 
interaction JT^]. Putting this together with Eq.(8) gives 

(AM) b5 = a s (2^ 5 )(mJm b ) - 3 [(AM) cc Va s (m c )]°- 65 [(AM) &5 /a s (m & )] - 35 . (10) 

It is interesting to note that for any reasonable choices of m c , and Aqcd, the quantity 
a s (2//b g )/a s (m c )°' 65 a s (m;,) ' 35 has a numerical value which is within about 3% of 1. (In a 
similar fashion, both QCD and relativistic corrections to Eq. (9) approximately cancel when 
fed into Eq.(10).) Moreover, (m&/m c ) ' 3 is always within a few percent of 0.7. Therefore, 
we can write 

(AM) bg ^ (0.7)(AM)«f(AM)f. (11) 
For the ground state, this reads 

Mb* - M Bc ~ (0.7) (Mj/^ - M„J - 65 (M T - M,J - 35 . (12) 
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Given the sources of uncertainty enumerated above, we expect this result to hold at about 
the 10% level. Only three of the six mesons appearing in Eq.(12) have been found experi- 
mentally thus far — namely, the J/ip, 1] C and T. Their masses are known quite accurately 
||16|| . There is some hope that the remaining three mesons will be detected in the near 
future at either the Fermilab Tevatron or LEP, allowing a test of the above result. 
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Potential 


System 


\R(0)\ 2 






IS 


2S 


3S 


4S 


5S 


6S 


Martin 


cc 


0.979 


0.545 


0.390 


0.309 


0.257 


0.222 




be 


1.720 


0.957 


0.685 


0.542 


0.452 


0.390 




bb 


4.423 


2.461 


1.763 


1.394 


1.164 


1.004 


logarithmic 


cc 


0.796 


0.406 


0.277 


0.211 


0.172 


0.145 




be 


1.508 


0.770 


0.524 


0.401 


0.325 


0.275 




bb 


4.706 


2.401 


1.636 


1.250 


1.015 


0.857 


Cornell 


cc 


1.458 


0.930 


0.793 


0.725 


0.683 


0.654 




be 


3.191 


1.769 


1.449 


1.297 


1.205 


1.141 




bb 


14.06 


5.681 


4.275 


3.672 


3.322 


3.088 


B uchmiiller- Tye 


cc 


0.794 


0.517 


0.441 


0.404 


0.381 


0.365 




be 


1.603 


0.953 


0.785 


0.705 


0.658 


0.625 




bb 


6.253 


3.086 


2.356 


2.032 


1.845 


1.721 


Lichtenberg- Wills 


cc 


1.121 


0.693 


0.563 


0.496 


0.453 


0.423 




be 


2.128 


1.231 


0.975 


0.846 


0.766 


0.711 




bb 


6.662 


3.370 


2.535 


2.139 


1.902 


1.740 



Table 1: Numerical values of the radial wave function at the origin squared, |-R(0)| 2 = 
|*(0)j 2 /47T, for the first six S-wave states of heavy quarkonium systems in various potential 
models. The parameters used in the potentials are discussed in the text. 
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